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Abstract. We present a summary of our previous investigations of the phys- 
ical and chemical structure of the Horsehead nebula, and discuss how these 
studies led to advances on the understanding of the impact of FUV radiation 
on the structure of dense interstellar clouds. Specific molecular tracers can be 
used to isolate different environments, that are more sensitive to changes in the 
FUV radiation or density than the classical tracers of molecular gas : the CO 
isotopologues or the dust (sub) millimeter continuum emission. They include the 
HCO or CCH radicals for the FUV illuminated interfaces, or the molecular ions 
H 13 CO+, DCO+ and other deuterated species (DNC, DCN) for the cold dense 
core. We discuss future prospects in the context of Herschel and ALMA. 



1. Introduction 

In this symposium, we have learned how the progresses in detector technolo- 
gies, with the development of Schottky mixers, InSb hot electron bolometers 
and finally SIS mixers led to spectacular discoveries on the chemical diversity 
and physical properties of the interstellar medium. The role of the interstellar 
medium in the life cycle of baryonic matter has been gradually uncovered in 
the last century. Following the discovery of the metallic lines in the spectra 
of reddened stars, the first molecules were recognized in the late 1930's. As 
these spectral features were extremely narrow, with velocity dispersions of a few 
kms -1 at most, the perspectives of radioastronomy and heterodyne spectroscopy 
for studying the interstellar medium, and especially its coldest and densest re- 
gions was quickly recognized. Molecular lines offer a wealth of diagnostics of the 
physical conditions, the chemical composition and the gas dynamics from the 
regions where they are formed. Molecular clouds themselves play a special role 
in the life cycle of interstellar matter, as the birth places of the new generation 
of stars but also as the main sites for interaction of stars with their environ- 
ment. Stars affect their environment through their radiative (through their X 
and UV radiation), mechanical (winds, shocks) and chemical (nucleosynthesis) 
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feedback. The understanding of the properties of molecular clouds, their struc- 
ture and their relation to the formation of stars has triggered the development of 
millimeter and submillimeter astronomy. Measuring the cooling radiation from 
molecular clouds is a key issue that determines our ability to detect the inter- 
stellar matter and understand how it can collapse and form stars. T. Phillips 
and his co-workers have brought a decisive contribution to this field, by mak- 
ing pioneering detections of the submillimeter lines of CO, atomic ca rbon and 
hydrid es. His prediction of the spectrum of a typical interstellar cloud [Phillips I 
(|l987l ) has been vastly cited and reproduced to illustrate the scientific discovery 
potential of submillimeter astronomy. 

The neutral interstellar gas cools through line emission, the most intense 
lines being the fine structure lines of ionized carbon, atomic oxygen and car- 
bon, together with the rotational lines of CO and water. Molecular hydrogen 
(H2) itself can be an important coolant, especially in dense photon-dissociation 
regions (PDRs) and shocks. Averaged over the whole interstellar medium of a 
galaxy, this line cooling radiation from PDRs makes a significant contribution to 
the total line emission at large scales dP hj Ums fc Huggins1ll98ll : iPhillips Ill987l : 
IWolfire et aT~l 120081 : ICubik et alH l2008h . iKeene et al. I (| 19851 ) have pioneered 
this field by mapping the atomic carbon line across the edge of two well known 
PDRs, S 140 and M 17. They showed that the atomic carbon line emission is 
very extended, and does not show the expected spatial pattern for a uniform 
density, plane parallel cloud. They discussed how the knowledge of the source ge- 
ometry and spatial structure are particularly important for understanding their 
line emission. 



2. Physical structure of the Horsehead nebula 



2.1. Summary of physical properties 

The Horsehead nebula is one of the most well known example of a molecular 
cloud. It appears as a dark patch against the bright emission from the HII region 
IC 434 at visible wavelengths, and as a bright line and continuum emission source 
at IR and (sub) millimeter wavelengths. Following the mapping of the CO(3- 
2) line emission with the CSO0, we have started investigating the physical and 
chemical structure of the Horsehead nebula. The source was chosen as it is 
a good template of photon dissociation region (PDR), where molecular gas is 
heated and photodissociated by Far UV (FUV; 6eV < hv < 13.6eV) radiation 
from nearby massive stars. The Horsehead nebula is particularly interesting as 
its geometry is simple. The illuminated edge is viewed very close to edge-on and 
the illuminating star (a Ori) lies at the same distance. The properties of the 
Horsehead nebula are summarized in Table [TJ 

Figures 1 and 2 present the li ne integrated in t ensiti es of the interstellar 



mole c ules we have observed so far (jTeyssier et al. 
20071 ): iGoicoechea et al. I (120061 . l2009h : iGerin et al. 



20041) : iPetv et al. I f|2005l . 



excited vibrational line of H2 at 2.1/im f rom lHabart et al. 
ISOCAM 7.7/xm map from lAbergel et al. I f|2003l ). 



2009)), t ogethe r with the 
" f|2005l ). and the 



1 see the CSO www page at www.submm.caltech.edu/ 
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Figure 1. Display of the maps obtained with the IRAM 30m and EfTelsberg 
(for NH3) telescopes towards the Horsehead nebula. All line maps show the 
integrated intensity in K.kms -1 between 9.0 and 12.0 kms -1 . The levels and 
color scale are shown on the right side of each map. The telescope beam is 
shown as a grey circle in the lower left corner. All maps have been rotated by 
14° counter-clock wise in order to display the PDR edge vertically. This edge 
is shown as the red vertical line. The white lines are drawn to help localize 
the main structures. 



2.2. Electron abundance 

The electron abundance ([e~] = n e /riH) plays a fundamental role in the chem- 
istry and dynamics of interstellar gas. It can not be directly measured in FUV 
shielded gas but can be derived from abundances and abundance ratios of key 
molecular species, which are sensitive to the presence of electrons. Stable molec- 
ular ions such as DCO + and HCO + have been traditionally used to estimate 
the electron abundance in cold cores because (i) they are abundant and eas- 
ily observable (ii) dissociative recombination is their main destruction route, 
and thus their abu ndances ar e roug h ly inversely proportional to the electron 
abundance (e.g. iGueljn et al I (| 19821 ): IWootten et al. I (| 19821 ): iMaret fc Bergin " 



(2007h: lHezareh et al. I f|2008l )). On the other hand, the presence of reactive ions 
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Table 1. Summary of the physical properties of the Horsehead nebula 



Name 

RA, DEC (2000) 

Distance from Earth 

Illuminating Star 

FUV Radiation field 

Depth along the line of sight 

Inclination of the line of sight 

Density profile 



Kinetic temperature 
Thermal pressure 



Barnard 33 (Horsehead nebula) 
05h40m54.27s, -02°28 / 00" 

- 400 pc (1" = 0.002 pc) 

aOri (09.5V) at 0.5° (3.5 pc), PA 76° 
60 (Draine unit) 

- 0.1 pc 
< 6° 

Steep gradient 

n ~ r~ 3 from ~ 10 5 cm -3 to 10 3 cm -3 in 10" 
and lower density halo 5 — 10 x 10 3 cm -3 
from Tk ~ 15 K in the cold core 
up to T K ~ 100K at the edge 

- 4 x 10 6 Kern" 3 



(species such as HOC + or CO + that react rapidly with H2) indica tes high elec- 
tron a b undance region s , e.g.. the UV irradiated cloud surfaces (e.g. JSmith et al."" * 



(120021 ); iFuente et al. I (|2003l )). The electron abundance is predicted to steeply 
decrease in a PDR from ~ 10 -4 at the cloud edge where the FUV radiation field 
is intense, down to < 10 -7 in the cold and well shielded gas. The magnitude and 
shape of this gradient is sensitive to the elemental abundances, especially sulfur, 
to the cosmic ray ionization rate, and to the efficiency of rec ombination reactions 
and electron attachment processes with PAHs an d grains dFlower et al. 1 120071 : 
IWakelam fc Herbst ll2008l : iGoicoechea et al. Il2009h . 

By co mbining pointed observat ions and maps of H 13 CO + , DCO + , HOC + 
and CO + , IGoicoechea et al. 1 ((2009) have investigated the electron abundance 
gradient across the edge of the Horsehead nebula. They showed that a con- 
sistent fit of both the PDR edge and the adjacent cold core can be obtained 
with a cosmic ray ionization rate of ( = 5 ± 3 x 10 -17 s _1 , assuming a low 
metal abundance of [M + ] = 10 -9 . The electron abundance in the Horsehead 
edge follows a steep gradient, with a scale length of ^0.05 pc (or ~25"), from 
[e~]~10 -4 times ^10 -9 in the core. Low metal abundances are required in the 
core, namely [M + ] < 4 x 10 -9 (with respect to H nuclei), to enable a significant 
deuterium fractionation of HCO + . This implies a very significant depletion of 
metals onto dust grains as the solar abundance of metals is [M + ] ~ 8.5 x 10 -5 ; 
( Anders fc Grevesse Ill989l ). The inclusion of PAHs in the chemical model mod- 



ifies the electron abundance gradient and decreases the metal depletion required 
in pure gas-phase models as a good agreement is obtained with metal abundance 
[M+] = 3 ± 1 x 10" 6 for a standard PAH abundance of 10" 7 . Indeed, PAH~ 
anions acquire large abundances, being more efficient in neutralizing atomic ions 
than electrons. PAHs have a significant impact on the negative charge budget 
of cold gas if their abundance is significant ([PAH]>10 -8 ). 
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Figure 2. Display of the images obtained with the IRAM PdBI and NRAO 
VLA (for NH3 ) towa rds the Horsehead nebula. The H2 2.1 /xm image is from 
lHabart et all (|2QQ5h and the ISOCAM 7.7 /xm image from lAbergel et al. I 
(|2003f ). Short spacings from the IRAM-30m (or Effelsberg 100m for NH 3 ) 
have been combined with the interferometer data to restore images with all 
spatial scales. The levels and color scale are shown on the right side of each 
map. The telescope beam is shown as a grey circle in the lower left corner. 
All maps have been rotated by 14° counter-clock wise in order to display the 
PDR edge vertically. This edge is shown as the red vertical line. The white 
lines are drawn to help localize the main structures. 



3. Chemical gradients 



3.1. Determination of chemical abundances 

Molecular line emission does not directly probe the abundances of a given 
species, even when several lines of the same species are detected. A minimum 
treatment of the molecular excitation and line radiative transfer is often needed 
to distinguish "excitation gradients" from true "chemical variations" . Our strat- 
egy is always to observe several rotational lines of the same species (sometimes 
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with line maps, sometimes using high signal to noise ratio spectra at represen- 
tative positions). In addition, we also try to observe molecular isotopologues, 
if possible, when the main line opacities can complicate the interpretation and 
lead to large errors on the derived abundance. Molecular excitation in PDRs is 
always far from LTE and different collisional and radiative effects compete at 
rates that depend on the physical conditions and on the molecule under study. 
In the last years, we have always tried to consistently link the observed line 
intensities with the abundances predicted by our chemical models by detailed 
non-LTE excitation and ra diative transfer modeling adapted to the Horsehead 
geometry (see Appendix in iGoicoechea et aTH (|2006l )). 

Table El and [3] list the molecular abundances for the PDR region and the 
dense core respectively, derived following the method described above. 



3.2. Observed chemical abundances 

It is obvious, when comparing the maps of various molecules shown in Fig. 1 and 
2 that the molecular abundances are not uniform across the Horsehead nebula. 
Molecular species like HCO and CCH closely delineate the western illuminated 
edge of the nebula. On the contrary, H 13 CO + and the deuterated species DCO + 
and DNC are mostly located in the dense cold core. Finally the dust continuum 
emission, as well as C 18 trace the bulk of the core, but do not show a sharp peak 
associated with the dense core. These maps have been discussed in detailed, and 
molecular abundances in our previous papers. 

In order to demonstrate quantitatively how these various species are spa- 
tially related, we show in Figure 3 a selection of scatter diagrams. The 1.2mm 
dust continuum and C 18 emissions are well correlated, while the correlation 
is less strong for H 13 CO + . The correlation between HCO and CCH is good 
as well, as both species trace the UV illumi nated matter in t he PDR- C4H and 



CCH are equally well correlated, as shown bv lTevssier et al. I (120041 ) : IPetv et al 



(120051). The recent confirm ation of the detection of CqH in the Horsehead edge 



by lAgundez et al. I ( 20081 ) confirms the high abundance of unsaturated hydro- 



carbons and radicals in this source. By contrast, the H 13 CO + (1-0) emission 
is tracing the dense cold core. It is loosely correlated with the dust continuum 
emission. The correlation of H 13 CO + with NH3 is quite good. 

The situation can be more complicated for other molecules. For instance, 
the C-C3H2 maps shows two emission regions, the PDR and the cold dense 
core. The presence of abundant C-C3H2 in the PDR is confirmed by the good 
correlation with HCO. The additional emission source shows up in the same 
plot as an upper envelope of points with strong C-C3H2 and relatively low HCO 
emission. This is confirmed in the C-C3H2 vs H 13 CO + scatter plot, that shows 
two clouds of points with two different slopes. Although most of the NH3(1,1) 
emission is produced in the dense and shielded region, some faint emission is 
detected in the PDR. This behavior shows up in the NH3 vs HCO and NH3 
vs H 13 CO + scatter plots as clouds of points deviating from the general trends. 
The ammonia data will be further discussed and analyzed in a forthcoming paper 
(Pety et al. 2009, in prep.). 
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Table 2. Summary of measured molecular column densities and abundances 
in the PDR. Offsets refer to RA(2000) = 05:40:54.27, DEC(2000) = -02:28:00 



Species 


HPBW 

[arcsec] 


Column density 
[cm -2 ] 


Abundance 

ATf V\ 
1M (A J 

N(H)+2N(H 2 ) 


Reference 


Offsets 
[arcsec] 


T I 


lz 


1.9 x ID 


U.o 


3,4,0,0 


(-5,0) 


O (J 


0.0 x 4.3 


4.U ± U.o x 1U 


1.9 x 1U 


I 


( A\ 

(-6,-4) 


r 1^ ITT 


7.Z X O.U 


o.U ± U.o X 1U 


1 a w 1 n — 8 
1.4 X 1U 


z 


( G. A\ 

(-6,-4) 


c-^3ri2 


1 v 4 7 


9441 n v 1 n 13 


1 1 v 1 n -9 

± . J. A ±U 


9 




c-C 3 H 2 


29 


9.3 ±0.2 x 10 12 


6.5 x 10- 10 


1 


(-10,0) 


I-C3H2 


27 


< 3.3 x 10 11 


< 4.6 x 10- 11 


1,7 


(-10,0) 


c-C 3 H 


28 


3.9 ±0.5 x 10 12 


2.7 x 10- 10 


1 


(-10,0) 


I-C3H 


28 


2.1 ±0.7 x 10 12 


1.4 x 10- 10 


1 


(-10,0 


C 4 H 


6.1 x 4.7 


4.0 ± 1.0 x 10 13 


1.0 x 10- 9 


2 


(-6,-4) 


C 6 H 


28 


9± 1 x 10 11 


2.2 x 10- 11 


8 


(-6,4) 


cs 


10 


8.1 ± 1.0 x 10 13 


2.0 x 10- 9 


3 


(+4,0) 


c 34 s 


16 


3.7 ±0.5 x 10 12 


9.2 x 10- 11 


3 


(+4,0) 


HCS+ 


29 


6.8 ±0.5 x 10 11 


1.7 x 10- 11 


3 


(+4,0) 


H 13 CO+ 


14 


5.8 ± 1 x 10 11 


1.5 x 10- 11 


5, 6 


(-5,0) 


HOC+ 


27.5 


1.8 ±0.6 x 10 11 


4.0 x 10- 12 


6 


(-5,0) 


CO+ 


10.4 


< 1.6 x 10 10 


< 2.7 x 10- 13 


6 


(-5,0) 


HCO 


14 


3.2 ±0.6 x 10 13 


8.4 x 10- 10 


5 


(-5,0) 



1 : iTeyssier et al. I (|2004l ); 2 : iPetv et~al71 (|2005f ). 3 : iGoicoechea et al I ([2j 
iPetv et al. I (|2007l), 5 : iGerin et al I (|2009l), 6 : IGoicoechea et al. I (|2009h 
iTevssier et aFI (j2QQ5h . 8 : lAgundez et"aFI (|2008l ) 



3.3. Chemical models of the Horsehead nebula 



Photo-dissociation models have been used to compute the physical and chemical 
structure of geometrically simple interstellar clouds, exposed to FUV radiations. 
Because of its well known geometry and FUV radiation field, the Horsehead neb- 
ula is a good template source for comparing model predictions with actual data. 
These comparisons have been extremely interesting in identifying the strong 
points of the present generation of models, and their remaining weaknesses, and 
in triggering new improvements. Figure 4 shows PDR model pre dictions adapted 
to th e Horsehead nebula, made with the Meudon PDR model dLe Petit et al 



20061 : IGoicoechea fc Le Bourlotl 120071 : iGonzalez Garcia et al 1120081 ). Predic- 



tions are s hown for the CO/C/C+ transition (jHabart ; et al. 1 120051). the sulfur 



chem istry (IGoicoechea et al. 1120061 ) . the electron abu ndance ( Goicoechea et al. 



l2009h. the deuter ium fr actionation flPetv et al. 1 120071). the HCO and CCH radi- 
cals flPetv et al. 1120051 : IGerin et al. T E)09) as well as for the main tracers of the 
nitrogen chemistry, CN, HCN, HNC and NH3. While the physical structure, and 
the transition from a FUV dominated chemistry in the PDR, to a cosmic ray 
dominated chemistry in the dense core, is well reproduced by the models, signif- 
icant discrepancies remain. The HCO, CCH and C4H radi cals are muc h more 
abu ndant in the PDR t han steady state model predictions. IPetv et aTH (|2005l ) 
and I Gerin et al. I ( 20091 ) have suggested directions for improvements. They in- 
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Table 3. Summary of measured molecular column densities and abun- 
dances towards the cold dense core. Offsets refer to RA(2000) = 05:40:54.27, 
DEC(2000) = -02:28:00 



Species 


HPBW 

[arcsec] 


Column density 
[cm -2 ] 


Abundance 

N(X) 
N(H)+2N(H 2 ) 


Reference 


Offsets 
[arcsec] 


H 2 


12 


2.9 x 10 22 


0.5 


6 


(+20,±22) 


CS 


10 


1.2 ±1.0 x 10 14 


2.0 x 10~ 9 


3 


(±21,±15) 


c 34 s 


16 


5.3 ±0.5 x 10 12 


9.1 x 10- 11 


3 


(±21,±15) 


HCS+ 


29 


6.8 ±0.5 x 10 11 


1.2 x 10" 11 


3 


(±21,±15) 


H 13 CO+ 


6.7 x 4.4 


6.5 ±2.5 x 10 12 


6.5 x lO" 11 


4,6 


(±20,±22) 


DCO+ 


12 


7.5 ±2.5 x 10 12 


8.0 x lO" 11 


4,6 


(±20,±22) 


HCO 


14 


4.6 ± 1.0 x 10 12 


8.0 x lO" 11 


5 


(±20,±22) 



1 : iTevssier et aTl (|2004 ): 2 : iPetv et all (|2005f ). 3 : iGoicoechea et al I ([20061 ), 4 : 
iPetv et al. I (|2007f ). 5 : iGerin et al.l (|2009h . 6 : IGoicoechea et al. I (|2009h 



elude i) the need to include more neutral-neutral reactions in the chemical net- 
works, that may have a significant impact on the chemistry (e.g. the O + CH2 
reaction in the HCO case), and ii) the need for considering the coupling of gas 
phase and solid phase chemistry as the photodesorption of molecular ices and/or 
the erosion of carbon grains may release significant amounts of radicals in the 
gas phase, that can modify current model predictions. 



4. Future prospects with Herschel, NOEMA, and ALMA 

Because of the steep gradients in the physical parameters and chemical abun- 
dances in PDRs such as the Horsehead nebula, an accurate comparison of ob- 
servations with models cannot be performed without information on the spatial 
variation of the molecular line intensities. Due to the small spatial scales in- 
volved, a combination of (sub) millimeter interferometer maps with short spac- 
ings from single dishes appears to be the best method for obtaining this spa- 
tial information at the required physical resolution. However, many interesting 
species show relatively weak lines, which are bearly detectable with the IRAM 
PdBI in a strong and nearby source like the Horsehead nebula (for instance the 
C4H lines). To access a more complete sample of molecular tracers, a more 
sensitive instrument is needed. With its larger number of antennas (12 vs 6), 
and wide instantaneous bandwidth (16GHz), the NOEMA project (NOrthern 
Extended Millimeter Array) for the IRAM Plateau de Bure Interferometer, will 
be able to produce high quality images of a much larger sample of molecules and 
sources, enabling the possibility to study the FUV chemistry in a wide diversity 
of physical conditions. These data will complement the information on the dust 
spectral energy distribution and gas thermal balance, which will be provided at 
unprecedented sensitivity and spatial distribution by ALMA. This will allow to 
probe for the first time the steep gradients in temperature and densities at PDR 
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Figure 3. Scatter diagrams of selected couples of images of the Horsehead 
nebula. All images have been resampled to the same spatial grid. Top line : 
The 1.2mm dust continuum versus C 18 0(2-l) and H 13 CO + (1-0). Left column, 
from top to bottom, scatter diagrams of CCH, C-C3H2 and NH3 with HCO. 
Right column : scatter diagrams of HCO, c-C 3 H 2 and NH 3 with H 13 CO + (l- 
0). 



edges, as well as the gradients in elemental abundances induced by freeze-out 
on grains. 
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Figure 4. PDR model adapted to the physical conditions in the Horse- 
head (x=60 and a density gradient see Goicoechea et al. 2009 for details). 
The uppermost panel shows the resulting H2 and H density profiles, to- 
gether with the expected gas temperature gradient. The next panels show 
the predicted abundances (with respect to to H nuclei) of the species we 
have observed and analyzed in detail in the past years: The CO/C/C + tran- 
sition (Habart et al. 2005) , The low sulfur depletion from CS observations 
(Goicoechea et al. 2006); The electron abundance, cosmic-rays ionization-rate 
and metal depletion (Goicoechea et al. 2009); The high HCO abundance and 
the HCO/H 13 CO + ratio as tracer un FUV-induced chemistry (Gerin et al. 
2009); the chemical deuteration of different molecules towards the cold and 
FUV-shielded regions (Pety et al. 2007); the small hydrocarbons emission in 
the PDR (Pety et al. 2005); the CN/HCN abundance ratio and the ammonia 
emission in the region (Pety et al. 2009 in prep.). 



Furthermore, hydrides play a key role in most chemical networks, as they 
are at the root of the formation pathways of most interstellar molecules (see the 
discussion by Lis et al. in this proceeding), but are however not accessible from 
the ground due to the large opacity of the Earth atmosphere at THz frequen- 
cies. With the full opening of the far infrared and submillimeter wavelength 
domain by the Herschel satellite, the spatial distribution of important hydrides 
such as water vapor H 2 0, OH, CH, CH+, NH, HF and many others including 
the isotopologues will be investigated. Although lacking spatial resolution, and 
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sensitivity, the pioneering detections of water by the ISO, SWAS and ODIN 
satellite have dem onstrated the added valu e of detecting hydrid e s for improving; 
the PDR models ([C ernicha ro et al.~1 120061 : iBergin et al. 1 120031 : iPersson et al. 



20091 ; iGonzalez Garcia et al. 1120081 ). The launch of the Herschel Space Obser 



vatory in May 2009 has concretized many decades of efforts for proposing a 
spatial submillimeter and far infrared telescope with direct detection detectors 
and heterodyne receivers. T. Phillips has been among the strongest advocates 
for such a space mission. The Herschel community warmly thanks him for this 
major achievement. 

Acknowledgments. M.G. and J. P. are indebted to T. Phillips for welcom- 
ing them in his group at CalTech. We enjoyed working in the CSO group, the 
stimulating atmosphere and Tom's insightful and enthusiastic comments on new 
ideas. The CSO is a wonderful telescope, where we learned many subtleties of 
performing heterodyne observations, from the behavior of local oscillators and 
mixers, to the pointing of the antenna ... 
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